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*Q | The various experimental data at AGS, SPS and RHIC energies on hadron particle yields for 

central heavy ion collisions are investigated by employing a generalized statistical density operator, 
that allows for a well-defined anisotropic local momentum distribution for each particle species, 
specified by a common streaming velocity parameter. The individual particle ratios are rather 
insensitive to a change in this new intensive parameter. This leads to the conclusion that the 
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reproduction of particle ratios by a statistical treatment does not imply the existence of a fully 
isotropic local momentum distribution at hadrochemical freeze-out, i.e. a state of almost complete 
thermal equilibrium. 
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I. INTRODUCTION AND MOTIVATION 



A principal question of the physics of (ultra-)relativistic heavy ion collisions is whether 
the strong interacting hadronic, or initially even deconfined partonic matter achieves an 
intermediate or final reaction stage of almost thermal and chemical equilibrium. One direct 
way to check for chemical equilibration is to consider and analyze experimental hadronic 
abundances data by fixing macroscopic intensive quantities such as temperature and chem- 
ical potentials within a statistical model of an equilibrated hadron gas The 
seminal idea of employing such a statistical description to hadron abundances goes back to 
Hagedorn for describing elementary reactions and to Hagedorn and Rafelski for the situation 
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of relativistic heavy ion co 
of Braun-Munzinger et al 



lisions [6j]. By now, numerous results, initiated with the work 
]], imply that a high degree of hadrochemical equilibration is 
achieved at a so called stage of chemical freeze-out, occurring in the central reactions of the 
heavy ion experiments carried out at the AGS at Brookhaven, at the SPS at CERN and at 
the RHIC at Brookhaven. 

The analyses are based upon rather restrictive assumptions (for two reviews having a 
critical look at the statistical model see [lUsl): The standard statistical model - often dubbed 
'thermal' model - assumes a complete chemical equilibrium (see jsj and references therein for 
allowing chemical non-equilibrium) as well as thermal equilibrium at hadrochemical freeze- 
out. The latter assumption implies the existence of isotropic momentum distributions locally 
at the individual space-time points within the fireball. The remarkable success of such 
'thermal' models popularized the belief that the applicability of such approaches also proves 
for having established a state of nearly complete thermal equilibrium, i.e. a state that has 
(locally) isotropic distributions in momentum space. In principle, however, a statistical 
description is more general. The statistical operator describes our knowledge about the 
system rather than the system itself, making use of incomplete information 10j. Still, to 
best of our knowledge, so far there has not been any investigation, concerning particle 
abundances, that starts from possible and appropriate statistical operators which do not 
simply resemble those of assumed (local) thermal equilibrium. This states the motivation for 
the present study, where we follow some general ideas of constructing statistical distributions 
based on information theory . Our picture is completely different from those that assume 
local isotropic momentum distributions supplemented by some hydrodynamical expansion 
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profile in space-time jll|. 

A possibility to probe the degree of local equilibration is to study the various stages of 



the reaction using microscopic transport models [12, 13, 141 (see also for lower bombarding 

nnn 

energies the work [la, [la, |17|). As ^ turns out, when testing 'thermalization', present 
transport models do not support the idea that a thermal system has been achieved at 
intermediate or even at later stages of the reaction, but show incomplete local momentum 
equilibration within the center of mass 

To meet this insight we propose in the following to relax the constraint of assuming lo- 
cally fully isotropic, statistical distributions for the comparison to hadron abundances, i.e. 
of assuming complete kinetic equilibrium. Certainly, at the onset of a heavy ion reaction, 
the local momentum distributions of the various particles are strongly non-isotropic. More- 
over, the various measured rapidity spectra point towards a global picture of the fireball 
which shows still a considerable longitudinal momentum excess as the spectra are much 
broader than those of a pure thermal and static source. Part of this broadening could be 
attributed to the buildup of a collective hydrodynamical expansion, but the other part can 
in principle originate from locally anisotropic momentum distributions. When constructing 
a generalized statistical operator, we will follow this latter possibility that the system might 
not have developed fully isotropic momentum distributions locally in space-time at the point 
of chemical freeze-out. 



II. GENERALIZED STATISTICAL OPERATOR FOR LONGITUDINALLY 
DEFORMED CONFIGURATIONS IN MOMENTUM SPACE 

A concept of a statistical description which incorporates at least partially the initial 
longitudinal momentum excess has been proposed by Neise ^| some time ago in order 
to allow for a generalization of the nuclear equation of state at finite density and temperatures 
for momentum anisotropic phase space distributions. Originally, the idea of two interpene- 
trating cold fluids goes back to Lovas [20j . who considered both an ellipsoidally deformed 
Fermi sphere and two separated Fermi spheres by introducing an additional parameter v 
which effectively alters the momentum distribution. This latter parametrization has been 
considered for various studies concerning the effective in-medium interaction and relativistic 
potentials experienced by nucleons inside a non-equilibrium environment 21] . For relativistic 
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heavy ion collisions at GSI energies, i.e. for Au+Au-collisions at 1 AGeV, the relative 
velocity v re i has been analyzed within a microscopic transport model and was found to 
be ~ 0.6 c at the stage of maximum compression manifesting a full nonequilibrium 
situation at intermediate reaction times. 

The main idea in is to interpret the velocity parameter as an additional statistical 
Lagrange multiplier v by introducing a generalized grand canonical potential 

Q = E - TS - fiN - vO , (1) 

where O represents the (thermodynamic) expectation value of an observable O, which the 
system should assume within the grand canonical average as an additional constraint. With 
this statistical approach various selfconsistency relations are automatically fulfilled. If the 
observable O is given by a one-particle operator, O = J2i n iOi, minimization of (JTJ) for a gas 



of independent hadronic partic 



the generalized distribution [Tja, |19j, |2 



es with a given set of one-particle energy states results in 
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drii ' exp [/3(ej — jJt, — uO{\ ± 1 



depending on the quantum statistics of the considered hadronic particle. The forma 
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can be easily generalized to allow for self-consitent mean-field type interactions [18, 
To allow for a specified momentum anisotropy, one can choose an adequate one-particle 
operator O characterizing the momentum anisotropy by taking the difference between the 
mean value of the momenta with a positive and negative longitudinal (z) component ^| 
respectively, i.e. 

°[p\ = EM(p^^(-ft)] fl P +fl p = J2\Pz\ a t a P ■ ( 3 ) 

V P 

According to (J2J one ends up with the generalized momentum distribution 

Hp = ex V [P(e p -v\p z \-^)]±l ' (4) 

with the conjugated intensive quantity v = v. 

A more appropriate relativistic notation for the distribution can be obtained by the 
following more intuitive and explicit construction of the same density operator [3]. Suppose 
that the homogeneous medium consists of two subsystems in momentum space (dubbed 
'projectile' and 'target' component indicating their origin), which are specified by their 



collective flow velocity i/ p )^ = 7(1,0, 0,i>) and v^ T ^ = 7(1,0,0, —v) in the combined 
center of mass frame with 7 = (1 — v 2 ) -1 / 2 . It becomes obvious that the intensive parameter 
v corresponds to the velocity of each subsystem. In the special cm. frame only the positive 
z-component, respectively negative z-component, of the (assumed) single-particle momenta 
contribute to the projectile subsystem, respectively target subsystem, i.e. both subsystems 
are separated by the plane p z = 0. For the grand canonical potential one takes the sum of 
the two components, i.e. 

n = v m» P w _ ToS w _ + v (t> p (t) _ TqS vo _ ^nP , (5) 

from which the statistical density operator immediately results as 

P = - exp {-fa (v^PP + vM'PF) + Po^Nb} . (6) 

For a noninteracting resonance gas of independent hadronic particles one has for the total 
momentum 

f m = EEW«(+/-ft) • ( 7 ) 

h P 

A similar expression holds for the entropy and baryon number. For the grand canonical 
potential (0) the contributions for the entropy and the baryon number are entering simply 
additively, and, as 6(p z ) + 0(—p z ) — 1, the explicit decomposition is only relevant for the 
sum v^ u Pj) P ^ + v( T ^ u P( T \ Putting everything together, from © one ends up with the 
expression 

K = — r r i — ( 8 ) 

exp (/Vy(ep - v\p z \) - PoPoj ± 1 
for the individual momentum distributions. (jHJ) is completely equivalent to (j4j) : fi and hq 
represent the true inverse 'temperature' and (baryo-) chemical potential defined in the local 
rest frame of each component, whereas (3 and /1 in (J3J are effective parameters modified by 7. 
A straightforward consideration shows that the momentum distribution becomes elongated 
in the longitudinal direction by the Lorenz factor 7. As an illustration, fig. depicts a 
situation typical for central Pb+Pb collisions at SPS energies for the case of nucleons, pions 
and kaons with an intensive parameter chosen to be v — 0.4 c. The heavier the hadronic 
particle, the more pronounced the two peaks ±m h/ ~fv in the distribution of the 'two fluids' 
manifest . 

Summarizing, the velocity v should be thought as the intrinsic (Lagrange-)parameter 
for characterizing any potential non-equilibrium situation in local longitudinal momentum 
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space within the statistical description. We consider it as the variable for characterizing a 
possible intermediate situation of the momentum distributions at hadrocemical freeze-out 
in the local comoving rest frame of a moving spatial cell within the spatially expanding 
fireball. This cell is then considered as a characteristic representative of the whole system, 
like for example within a boost-invariant Bjorken geometry: A boost-invariant but locally 
momentum anisotropic distribution is given by the standard substitution 

= n h (p ± ,p z ) = n h (p u y) -»■ f h (p ± ,y;z,t) = n h (p ± ,y-r]) (9) 

of Lorentz-boosted comoving cells, and assuming the intrinsic parameters T, \lb and v only 
to depend on the eigentime r. 

The parameter 'jv is certainly different and considerably smaller than in the very initial 
situation of the heavy ion reaction due to the (unknown) amount of dissipation occurring in 
the (violent) initial stages. The mismatch compared to a fully momentum isotropic situation 
should relax continuously in time due to collisions among the various particles, although the 
system might never actually become locally completely isotropic Q,^^, as outlined in the 
introduction. In the following we thus will consider v as a free parameter. 

Finally it is worthwhile to point out the strict meaning of a thermodynamical equilibrium 
description compared to statistical descriptions in general [10]. A statistical operator, relying 
on the concept of incomplete information, assumes the existence of some average constraints 
(Of). The density operator then has the standard form p = exp(fi — J2% v iOi) [10]. Generally, 
the operator Oi does not need to commute with the Hamiltonian H. Being more restrictive, 
one has a stationary ensemble when taking only constants of motion for the Of, in other 
words, the Liouville equation for the density matrix 

f f ={H,p} = 0, (10) 

has to be fulfilled for a true thermodynamical equilibrium situation. This is not the case 
for the chosen momentum anisotropy operator O reflecting the fact that the momentum 
anisotropy defined might actually be a time dependent quantity. As outlined above, col- 
lisions among the particles in the (strongly) interacting matter are ultimately destroying 
the stationarity. A requirement for stationarity is fulfilled when considering spatially homo- 
geneous matter and employing approximate effective one-particle (and thus collision-less) 
Hamiltonians of mean-field type or a simple quasi-free Hamiltonian for describing the hadron 



resonance gas. On the other hand, the statistical assumption is made when postulating the 
existence and minimization of the generalized grand canonical potential (0) with the con- 
straint that locally a finite mismatch in momenta survives at hadro-chemical freeze-out. 



III. INFLUENCE OF THE VELOCITY PARAMETER ON PARTICLE RATIOS 



Following the bootstrap idea of incorporating the strong interactions among the con- 
stituents (see eg [s^)' we assume a 'noninteracting' hadron (Hagedorn-like) resonance gas 
for describing the excited hot hadronic matter inside the fireball [l|, . The chemical 

freeze-out of particles is assumed to happen for a fixed temperature T. For the sake of 
simplicity we do not intend to distinguish between a 47r-analysis or an analysis of particle 
ratios being measured only close to midrapidity: An approximate Bjorken geometry as given 
in (JHJ) would give the same answer for both analyses. 

We choose the description of an isospin symmetric hadronic gas and include baryonic and 
mesonic resonances of masses up to 2 GeV 

For the various particle number densities 
described by the generalized momentum distributions (JBI) one has 



Pi {Lf / rf3p exp(/3 o7 (^ - v\ Pz \) - A^o) ± 1 ' (H) 



with spin degeneracy gi, Po = and Ei = y p]_ + p 2 z + mf. The chemical potential of each 
hadronic particle is specified by its baryon number content Nb and its strangeness number 
content N$ via /ij ; o = Nb^b + Nsfis, assuming isospin symmetric hadronic matter (/xj = 0). 

Changing the velocity parameter, the various particle number densities do increase with 
increasing v to almost 2 / y for velocities close to c. The total intrinsic energy density, on the 
other hand, will scale with almost 2 / y 2 for such large velocities. 

We are now in position to check for the sensitivity in extracting the thermodynamical 
intensive parameters like the 'temperature' parameter (3q 1 and the chemical potentials when 
varying the velocity parameter v. It is worthwhile to first discuss two limiting situations: 
For v — the standard statistical model description results. For velocity parameters being 
sufficiently large, the two interpenetrating components fully separate in momentum space, 
so that the hadronic particle ratios will not be altered when employing the same temperature 
Pq 1 and chemical potentials as for v — 0. Each of the separate subsystems then describes the 
situation of a Lorentz boosted momentum isotropic fireball, i.e. the situation for v = 0. The 
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question remains for the true intermediate situations, where one cannot really distinguish 
among two separate components. As a further limiting case, employing a classical Maxwell- 
Boltzmann approximation, 

n p « exv(-(3 j(Ei-v\p z \)) exp(/3 /^, ) , 

number ratios of particle and their anti-particles are not altered by a change in the velocity 
parameter, being solely dependent on the exponential factor containing the chemical poten- 
tial. We therefore expect only some smaller kinematical sensitivity for particle ratios with 
unequal masses. 

For testing the potential influence of an anisotropic momentum distribution on the 
hadronic particle ratios the strategy is as follows: We first employ the known values for 
the temperature T and potential fiB,o of Q, 0, 0, 0] within the present description taking 
v — for describing the various situations of heavy ion experiments carried out at the 
AGS, at the SPS and at RHIC. As one particular example we show in fig. El characteristic 
particle ratios in comparison to the experimental results for central Pb+Pb collisions at 
SPS energies. We find_good agreement to the various hadronic particle ratios (in quality 
like those given in Q, JO, Q|) and we mention the well-known fact that feeding from all 
relevant hadrons and resonances to the abundances of 'stable' (i.e. detectable) particles is 
of crucial importance. With the temperature and the chemical potentials held fixed, we 
now change the velocity parameter in order to see the sensitivity on the particle ratios. In 
the following we concentrate on the p/p - , the p/ii + - and the K + /K~ -ratio. The first 
two ratios typically fix the temperature and the baryochemical potential. The last fixes the 
strange chemical potential, which, however, enters as a constraint as the system carries no 
overall net strangeness. In fig. El the actual (no n-) dependence of these exemplary particle 
ratios on varying the velocity parameter is depicted. As expected from the above estimates, 
the largest deviation, about 6%, one finds for the p/ir + ratio (for the situation of a Pb-Pb 
collisions at SPS). We have tested for all the various other particle ratios shown in fig. |21 
and find, if at all, only a marginal dependence on the velocity parameter v. p/ir + shows 
the 'strongest' sensitivity. The same situation is met when looking at ratios measured at 
AGS or at RHIC. Most of the other ratios, such as the shown p/p and K + /K~, only vary 
by about 1% or even less. In tab. 0the above three particle ratios as obtained within our 
investigation are summarized. There is no remarkable change of the particle ratios at any 
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of the considered collision energies. 

Hence, the quality of a possible detailed fit to the particle ratios including the velocity as 
a free parameter would not be changed: An extraction of v would be completely ambiguous. 



IV. SUMMARY AND CONCLUSIONS 

We have constructed a statistical description of a hadronic resonance gas allowing for 
a finite longitudinal momentum excess compared to commonly employed isotropic ('ther- 
mal') distributions. This is achieved via a velocity parameter acting as a common Lagrange 
multiplier for fixing a momentum anisotropy for the local momentum distributions . As 
microscopic transport models do not fully support the idea of a complete local momentum 
equilibration at any stage of a central heavy ion reaction, the presented approach relaxes 
this typical constraint of a standard statistical model. From a physics point of view one 
can argue that the amount of dissipation or collisions is simply not enough to achieve full 
momentum equilibrium locally in space. Going further with the argument, there exists al- 
ways the principal question when applying a statistical treatment for hadron production 
whether it represents a reflection of a true dynamical and stationary equilibrium state or, 
more generally, whether it is only a reflection of a state minimizing the content of informa- 
tion (i.e. maximizing the 'entropy') within a set of appropriate, global constraints [10], as 
also employed in this study. This is particularly obvious when applying such concepts to 
e + e - -, pp- or pp-collisions, as proposed seminally by Hagedorn and described recently by 



Becattini 2J| within a micro-canonical treatment. For such 'systems' one can not expect 
the occurrence of a true dynamical equilibrium state (fulfilling approximately ([10))) at any 
stage of the reaction. 

It turns out that the hadronic particle composition is more or less completely insensitive 
to varying the velocity parameter from zero to nearly speed of light. The extraction of this 
parameter can not be achieved by looking at particle number ratios. On the other hand, 
the extraction of the 'temperature' parameter as well as the (baryo-) chemical potential, 
however, are the same as with the standard approach. They are thus the significant and 
appropriate statistical quantities for dealing with the various number ratios characterizing 
the intermediate stage of the reaction at hadrochemical freeze-out. 

The popular conclusion that the many successes in describing particle ratios within a 
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'thermal' model give evidence for the existence of a state of thermal equilibrium at hadro- 
chemical freeze-out is, however, not true. Local isotropic momentum distributions reflect 
merely a simplifying assumption of the applied statistical operator. With this observation 
in mind, the standard 'extraction' of total particle number densities or energy densities, 
respectively, in absolute terms achieved inside the hadronic fireball at chemical freeze-out 
becomes also unsteady. Our study underlines the necessity of understanding relativistic 
heavy ion collision experiments by means of detailed microscopic transport models. 
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model for v = 


model for parameter v 
with largest deviation in % 


experiment 


AGS 


T = 120 MeV, [is = 540 MeV, ptg = 108 MeV 


P/p 
p/n + 
K+/K- 


0.000124 
1.0535 
5.5362 


0.04652 
5.9178 
-0.5315 


0.00045 ± 0.00005 
4.4 ±0.4 


SPS I 


T = 160 MeV, fi B = 170 MeV, fx s = 38 MeV 


P/p 
p/ir + 
K+/K- 


0.1195 
0.1863 
1.51 


0.0048 
6.1940 
-0.2680 


0.12 ±0.02 
0.18 ±0.03 
1.67 ±0.15 


SPS II 


T = 168 MeV, \i B = 266 MeV, /x 5 = 71.1 MeV 


p/p 
p/ir + 
K + /K~ 


0.0422 
0.3294 
2.0820 


0.0104 
4.7817 
-0.4319 


0.055 ±0.01 
1.85 ± 0.09 


RHIC 


T = 174 MeV, fi B = 46 MeV, ^ s = 13, 1 MeV 


P/P 
p/ir + 
K+/K- 


0.5845 
0.1129 
1.1282 


0.8568 
5.4299 
-0.1297 


0.60 ± 0.07 
1.136 ±0.070 



TABLE I: Characteristic particle ratios from various heavy collision experiments (from top to 
jottom): (1) Si + Au at the AGS, parameters from |lj; (2) S + Au at the SPS, parameters from 
3); (3) Pb + Pb at the SPS, parameters from _3]; (4) Pb + Pb at RHIC, parameters from Q. The 
maximal deviation (in percent) by changing the velocity parameter is indicated. 
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Pz/T " V 20 Pt/T 




v=0.4 

FIG. 1: Generalized statistical momentum space distribution with v = 0.4 for the nucleons, the 
pions and the K + . Temperature parameter (3q and chemical potentials are taken as those extracted 
for v = being typical for Pb+Pb collision at SPS energies Qj. 
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FIG. 2: Characteristic particle ratios for Pb+Pb collisions at SPS energies obtained with velocity 
parameter v = 0. Temperature parameter (3q and chemical potentials are taken the same as in 3] 
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FIG. 3: (Weak) Dependence of three typical particle ratios (p/p; p/tt + ; K + /K~) varying the 
velocity parameter from v = c to v = 0.9 c. Temperature parameter and chemical potentials are 
taken the same as those of fig. |21 with v = 0. 
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